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A generalized Anderson model for a magnetic ion in a harmonic potential is formulated. 
The model is investigated by the numerical renormalization group(NRG) method. In addition 
to the conventional s-wave screening, the model exhibits phonon assisted p-wave Kondo effect 
as well as Yu and Anderson type Kondo effect. It is shown that the s-wave Kondo and the 
Yu- Anderson Kondo belong to the same fixed point. At the boundary between the s-wave and 
p-wave Kondo regions line of fixed points of the two channel Kondo effect is identified. 

KEYWORDS: Kondo effect, impurity Anderson model, electron-phonon coupling, numerical renormaliza- 
tion group, two channel Kondo effect 



Filled skutterudite compounds RT4X12 (R= rare earth 
or alkaline earth element; T=Fe, Ru, Pt, or Os; X=P, As, 
Ge, or Sb) are characterized by their specific structures 
which involve a network of cages filled by guest ions. 
When the radius of the guest ion is smaller than the di- 
ameter of the cage, the guest ion vibrates with larger 
amplitude and smaller frequency than conventional lo- 
calized modes. Such anharmonic local vibrations are re- 
ferred to as rattling modes and may lead to novel phe- 
nomena of the coupled electron-phonon systems. Recent 
experiments on SmOs4Sbi2 show a large electronic spe- 
cific heat coefficient linear in T, which is robust against 
magnetic field. 1 -* One possible scenario of the unusual be- 
havior is the effect of strong electron-phonon coupling. 

When a magnetic ion couples with conduction elec- 
trons we expect Kondo effect due to the localized mo- 
ment. On the other hand, it has been shown by Yu and 
Anderson(YA) that when conduction electrons couple 
strongly with ionic vibrations a different type of Kondo 
effect is expected. 2 ' In this Letter we will study the in- 
terplay between the conventional Kondo effect and the 
Yu-Anderson type one. Concerning the effects of cou- 
pling between a magnetic ion and ionic vibrations, Hotta 
has studied effects of anharmonicity in the Holstein- 
Anderson model. 3 ) A lattice version of the Holstcin- 
Anderson model is also studied. 4 ) In this Letter, we will 
discuss effects of transverse vibrations on the Kondo ef- 
fect rather than the breathing type vibrations. 

Suppose electron orbitals of a magnetic ion are de- 
scribed by the £-spherical-wave functions. Then the hy- 
bridization with the conduction electrons is given by 
overlap integrals between the £-wave localized orbitals 
and plane waves. To include the effect of vibrations the 
origin of the localized orbitals is shifted by Q which is the 
coordinate of the ion position. We can expand the overlap 
integral with respect to Q. In the zeroth order the local- 
ized ^-spherical- waves hybridize with the ^-partial-waves 
of conduction electrons. In the first order of Q, they hy- 
bridize with the I + 1 and £ — 1 partial waves. In the 
present Letter we will concentrate on the simplest case 
of s-wave localized orbital. Then the total Hamiltonian 
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up to the order of Q is given by the sum of 
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where i=x,y,z represent three components of displace- 
ment vector as well as the three components of £=\ par- 
tial waves in real representation. In this paper phonons 
are assumed to be harmonic, and in Eq.(l) <Zj is an an- 
nihilation operator for the i-component phonon and uj is 
the frequency of the phonons. Annihilation and creation 
operators for the localized electrons are described by j a 
and fl, ef is the energy of the localized orbital and U 
the Coulomb interaction. Vq is the matrix element of the 
hybridization of the s-wave channel and concerning the 
p-wave channel hybridization takes place assisted with 
phonons with the matrix element V\. 
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Fig. 1. (Color online) (a) NRG scheme of the present model. 
The wavy line represents hopping processes accompanied with 
phonon creation or absorption, (b) The NRG scheme for the 
original Yu-Anderson model, (c) Low energy fixed point of the 
conventional s-chK. (d) Low energy fixed point of the p-chK. 



To investigate properties of the model we first simplify 
the model to include vibrations of only one direction, say 
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i—x. We use the NRG method to analyze the model in 
detail. 5 ' 6) In the NRG algorithm we assume a constant 
density of states p=\/2D with band width of 2D and 
discretize it in logarithmic energy scales controlled by 
A for both the s-wave and p-wave channels. Then the 
present problem is mapped to the set of Hamiltonians 
where the localized orbital is coupled with two Wilson 
chains corresponding to the s-wave and p-wave channels 
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where the quantities with tilde are the parameters 
modified by the logarithmic discretization. Figure. 1(a) 
shows the present system schematically. We assume the 
electron-hole symmetry and thus the average number of 
electrons is kept to be one per site. 

We first discuss several limiting cases. When the 
Coulomb interaction is strong electrons in the mag- 
netic ion behave as a magnetic moment. When Vo is 
stronger than V\, the magnetic moment is quenched 
mainly through the coupling with the s-channel conduc- 
tion electrons as temperature is lowered. This situation 
is a case of standard Kondo effect and will be referred 
to as s-channel Kondo (s-chK) regime in this Letter. On 
the other hand, when V\ is dominant over Vo, the mag- 
netic moment is screened mainly by the p-wave channel 
and this regime will be called p-channel Kondo (p-chK) 
scheme. 

When both the Coulomb interaction and the phonon 
assisted hybridization with the p-wave channel are not 
important the system may be treated as a renor- 
malized Fermi(RF) chain. An interesting case in the 
weak Coulomb interaction regime is realized when the 
electron-phonon coupling V\ is strong. It is instructive 
to consider the two-site problem which is described by 
the following Hamiltonian: 

# 2 -sitc = Vi [pLu + ( h - c -)] ( a + at ) + ^ at « ( 5 ) 

a 

This problem is solved easily by the canonical transfor- 
mation, 7 ) e lS with S = i\EA~ A L A o* + Bla B 0a){a - 
a') where Aq g and Bq„ are the annihilation operator for 
the antibonding and bonding orbitals, A§ a = -^(/<t — 
Poa) and B 0a = -^(/<r + Poa) -The coupling constant 
is given by A = —Vi/fko. The ground states are doubly 
degenerate 

\L)=B t ot Bl i \™c)\£) with \£)= e- 2A ( a ~ at )|0) (6) 
|R)=4 t 4;l vac >l r ) with |r)^e 2A ( Q - at )|0) (7) 

and the binding energy is Epv = — 4 1+ ^_ t . In the 
state |L), the bonding state is doubly occupied by up 
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Fig. 2. (Color online) (a) Entropy of the noninteraction Yu- 
Anderson model, (b) Temperature dependence of entropy for 
various value of phonon frequency u. In the panel (a) the red 
squares are the points where the entropy is larger than 1 . Ike log 2 
and the green dots are for 0.9kg log 2 < Si mp < 1.1kg log 2, and 
the blue triangles for Si mp < 0.1kg log 2 with the fixed points 
corresponding to Fig. 1(c). 

and down electrons forming a spin-singlet state and the 
ionic position is shifted to left, while in the state |R), 
the antibonding state is doubly occupied and the ionic 
position is shifted to right. 

Kondo effects discussed by YA represent screening pro- 
cess of the |L) and |R) states. In the original model of 
YA, the localized orbital is not introduced and only scat- 
tering processes of conduction electrons between the s- 
wave and p-wave channels are considered. Therefore the 
original model by YA may be expressed by Fig. 1(b). 
However in the non-interacting case, U=0, essential 
physics of the two models is the same. Therefore the 
model discussed in this Letter spans the Kondo effects 
of a magnetic impurity to the YA type Kondo effects 
(YAK) due to strong electron-phonon coupling. It should 
be mentioned that the same model as the present one is 
discussed by Dias da Silva and Dagotto in connection 
with phonon-assisted tunneling in molecular junctions. 8 ) 

In this Letter, we will show first that YAK are ac- 
tually realized in the present model by using the NRG 
and then determine the phase diagram in the U-V± plane 
which clarify the essential physics of the Kondo effects 
of a vibrating magnetic ion. 

Now we will show results of NRG calculations. Con- 
cerning the parameters of the NRG, we keep 15,000 
states with 50 phonon states and use the cutoff A=3.0. 

Let us start from the non-interacting case. Fig. 2(b) 
shows temperature dependence of entropy for the impu- 
rity part including entropy associated with the phonon 
degrees of freedom. For this example the coupling con- 
stants are kept fixed, Vb=Vi = 0.2 and the phonon fre- 
quency uj is varied. When uj gets smaller more phonons 
are excited, leading to effectively stronger electron- 
phonon interaction. For w=0.300 the entropy is monoton- 
ically released as temperature is lowered, which suggests 
that the system is in the RF chain regime. As u decreases 
one can see emergence of a plateau in the entropy curve 
at around ke log 2. This is due to formation of polaron 
doublet(PD) which correspond to |L) and |R). One can 
see that entropy release after the plateau for u;=0.262, 
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Fig. 3. (Color online) Effect of Coulomb interaction for the Yu- 
Anderson model, (a) Phase diagram of in U-T plane for U > U c , 
where U c = 0.00194 is the Coulomb interaction where a 2chK 
is observed, (b) For U < U c - (c) Temperature dependence of 
entropy for various U > U c . (d) For U < U c . In the panel (a) 
and (b) the purple triangles(downward) are the points where the 
entropy is 0.4kelog2 < Si mp < 0.6kBlog2, and the aqua dia- 
monds for Sj mp < O.lkB log 2 with the fixed points corresponding 
to Fig. 1(d). 



0.245, and 0.238 are typical examples of the YAK. To 
the best of the authors' knowledge this is the first NRG 
calculation for the spinful YAK. When we proceed to 
the four-site problem from the two-site problem one can 
see that matrix elements which connect the bonding and 
antibonding orbitals appear. Therefore, in the present 
model matrix elements between the |L) and |ft) naturally 
appears on this level and YAK are found in a reasonable 
range of parameters. On the contrary, when we use the 
original YA model, Fig. 1(b), and assume complete sym- 
metry between the s- and p-wave channels, there is no 
matrix elements for the mixing. Therefore higher order 
processes like two-phonon processes are required for the 
YAK to be realized. 9 ' 

When u> gets smaller, then the binding energy of the 
PD becomes larger and the Kondo temperature of the 
YAK becomes extremely small. Therefore YAK are not 
seen in the temperature range studied when u> < 0.235. 
The panel (a) summarizes the results on the entropy in 
U-T plane. One can estimate the Kondo temperature of 
YAK from the middle of the white region. 

To sec the effects of Coulomb interaction on the YAK, 
finite U is introduced for a fixed w=0.245. The results 
are shown in Fig. 3. The panels (b) and (d) show that 
when the Coulomb interaction is weak the behaviors of 
the YAK do not change. However, in a narrow range of U 
around U=0. 00194 a new plateau appears in the temper- 
ature dependence of the entropy. The new plateau region 
is shown by purple triangles (downward) in the panels 
(a) and (b). The fact that the entropy at the plateau is 
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Fig. 4. (Color online)The phase boundary between the s-chK 
fixed point and the p-chK fixed point. 



(l/2)kBlog2 suggests that the plateau is due to a two- 
channel Kondo effect(2chK). Actually our analysis of the 
energy spectrum at this critical U confirms that the new 
fixed point is the 2chK fixed point. 10 ' Through the crit- 
ical U c =0. 00194 the plateau region of ke log 2 at higher 
temperatures remains for U < 0.2 which is of the order of 
the binding energy of the PD estimated to be 0.24. We 
have looked susceptibility of the system and it is con- 
firmed that the plateau is due to the PD since there is 
no enhancement of the susceptibility in this parameter 
range. 

Dias da Silva and Dagotto have found a 2chK behav- 
iors for the set of parameters Vb=0.2, Vi =0.1685, [7=0.5, 
w=0.2. 8 ' We have confirmed that the 2chK behavior is 
seen for the set of parameters Vb=0.2, 14=0.1646, [7=0.5, 
and w=0.2. The small difference of V\ can be attributed 
to difference in details of numerical calculations. It is im- 
portant to figure out the relation between the 2chK fixed 
point in the intermediate coupling region and the one in 
the close vicinity of the YAK. 

Figure. 4 shows the separatrix on which the 2chK be- 
haviors are observed in the U-V\ plane for Vq= 0.2, oj =0.2 
fixed. When we consider the situation where typical s- 
chK is observed the energy spectrum for finite systems 
at the low energy fixed point shows oscillatory behav- 
iors consistent with the schematic Fig. 1(c). On the other 
hand when the system shows p-chK, oscillatory behav- 
iors of the energy spectrum is characterized by Fig. 1(d). 
We have found for any set of parameters the energy spec- 
trum shows one of these behaviors. The separatrix is the 
boundary between the two behaviors. 

The upper left corner of Fig. 4 corresponds to s-chK 
and the upper right corner to p-chK. Therefore it is clear 
that the boundary in the large U region represents the 
2chK whose origin is competition between the s-wave 
screening and the phonon assisted p-wave screening. The 
degrees of freedom to be screened is the spin of the lo- 
calized moment on the magnetic impurity. 

Now we decrease U from the region of the s-chK for 
Vi = 0.164 fixed, Fig. 5. For U = 1.2, we see a typical 
Kondo behaviors of the entropy. Around [7=0.5, in low 
but finite temperature ranges, we see reminiscent of 2chK 
behaviors which show that the finite temperature prop- 
erties are influenced by the unstable 2chK fixed point. For 
small [7 the RF chain behaviors are observed. 
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Fig. 5. (Color online)Temperature dependence of entropy for sev- 
eral parameters along the cut Vi =0.164 in Fig. 4. 

On the other hand along the cut of Vi=0.1666, Fig. 6, 
we first observe the s-chK behaviors for £7=1.5 with a 
shoulder corresponding to Icb log 2 due to the formation 
of a magnetic moment. Then at U =0.99 the first 2chK 
behaviors is observed. At U =0.8 we observe a p-chK be- 
havior influenced by the 2chK fixed point. As we ap- 
proach the lower boundary at U =0.25 we see typical 
behaviors of the 2chK again. For [7=0.0 a shoulder at 
ke log 2 is seen which is due to the PD. From the phase 
diagram we have established that the 2chK behaviors in 
the strong correlation regime is continuously connected 
to the 2chK in the close vicinity of YAK in the weak 
correlation regime. 

The latter 2chK is associated with the PD. From the 
point of view of electron states one component of PD, at 
the two-site level, is the bonding singlet while the other is 
the antibonding singlet as shown in eq.(6) and (7). These 
two degrees of freedom may be also described by -0ls = 
-^(flft +P T H)|vac) and ^ps = ^ 2 {f\p\ ~ /> T )|vac) 
which represent local singlet channel and pair singlet 
channel. For the noninteracting case of the two site prob- 
lem they have precisely the same weight resulting in the 
PD. When we attach the s-channel and p-channel chains, 
the two local singlets will be extended into chains to 
gain kinetic energy, leading eventually to the same fixed 
point as the s-chK. On the other hand when U is fi- 
nite the pair singlet channel becomes more relevant. At 
the critical point the two channels exactly balance and 
the 2chK takes place. Important point is that concerning 
the YAK, the degrees of freedom to be screened are the 
ion displacement described by \t) and |r) and electronic 
part defines the channels to screen. The two different 
channels are schematically shown by Fig. 1(c) and (d). 
At this point we would like to comment on the theory 
by Kusunose and Miyake 11 - 1 who proposed a possibility 
of the 2chK for the YAK. They considered that the spin 
degrees of freedom of conduction electrons played the 
role of screening channels. However, as the wave func- 
tions |L) and |R) show, electronic parts of the PD are 
spin-singlets and their discussion is not justified for the 
spinful YAK. 

In conclusion we have derived a generalized Anderson 
model for a vibrating magnetic ion coupled with con- 
duction electrons. We have studied the simplest case of 
s-wave localized orbital and ion vibrations in one direc- 
tion. The model shows the s-chK and p-chK as well as 
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Fig. 6. (Color online) Temperature dependence of entropy for sev- 
eral U along the cut V\ =0.1666. 

the YAK behaviors. The relation among them is clarified 
by the NRG calculations. It is shown that there is a 2chK 
unstable fixed points in the close vicinity of the YAK in 
weak correlation regime. The boundary is continuously 
connected with the 2chK behaviors in the strong corre- 
lation regime. 

The present results suggest that Kondo behaviors, in 
particular a large mass enhancement, may be realized in 
a wide range of parameters once the vibration frequency 
of a magnetic ion is small and the electron-phonon cou- 
pling is strong. This mechanism of mass enhancement 
takes place in the weak correlation regime and therefore 
is expected to be robust against magnetic fields. In this 
Letter we have treated only harmonic phonons. However, 
there is a growing consensus that rattling type anhar- 
monic phonons provide low energy vibration modes with 
effectively strong coupling. 12 ~ 14 ) It is an interesting fu- 
ture problem to investigate Kondo effect of a vibrating 
magnetic ion in strongly anharmonic potential. 
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